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Abstract 

Nucleon-nucleon (NN) cross sections are evaluated in neutron-rich matter 
| using a scaling model according to nucleon effective masses. It is found that 

o ' 

the in-medium NN cross sections are not only reduced but also have a dif- 
ferent isospin dependence compared with the free-space ones. Because of 
the neutron-proton effective mass splitting the difference between nn and pp 
scattering cross sections increases with the increasing isospin asymmetry of 

X 

$— i ' the medium. Within the transport model IBUU04, the in-medium NN cross 

sections are found to influence significantly the isospin transport in heavy- 
ion reactions. With the in-medium NN cross sections, a symmetry energy 
of E sym (p) f» 31.6(p/ po) om was found most acceptable compared with both 
the MSU isospin diffusion data and the presently acceptable neutron-skin 
thickness in 208 Pb. The isospin dependent part K asy (po) of isobaric nuclear 
incompressibility was further narrowed down to —500 ± 50 MeV. The pos- 
sibility of determining simultaneously the in-medium NN cross sections and 
the symmetry energy was also studied. The proton transverse flow, or even 
better the combined transverse flow of neutrons and protons, can be used as a 
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probe of the in-medium NN cross sections without much hindrance from the 
uncertainties of the symmetry energy. 
PACS numbers: 25.70.-z, 25.75.Ld., 24.10.Lx 
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I. INTRODUCTION 



The isospin dependence of in-medium nuclear effective interactions is important for un- 
derstanding not only novel properties of exotic nuclei near drip lines but also many in- 
teresting questions in astrophysics [1-3]. Especially, it determines the symmetry energy 
E sym (p) term in the equation of state (EOS) of isospin asymmetric nuclear matter. The 
density-dependent symmetry energy itself is still poorly known but very important for both 
nuclear physics and astrophysics. Heavy-ion reactions induced by neutron-rich nuclei pro- 
vide a unique opportunity to explore the isospin dependence of in-medium nuclear effective 
interactions, especially the symmetry energy, in a broad range of density. This is because the 
isospin degree of freedom plays an important role in heavy-ion collisions through both the 
nuclear EOS and the nucleon-nucleon (NN) scatterings [4,5]. In particular, the transport of 
isospin asymmetry between two colliding nuclei is expected to depend on both the symme- 
try potential and the in-medium NN cross sections. For instance, the drifting contribution 
to the isospin transport in a nearly equilibrium system is proportional to the product of 
the mean relaxation time r np and the isospin asymmetric force F np [6]. While the r np is 
inversely proportional to the neutron-proton (np) scattering cross section a np [6] , the F np is 
directly related to the gradient of the symmetry potential. On the other hand, the collisional 
contribution to the isospin transport in non-equilibrium system is generally expected to be 
proportional to the np scattering cross section. Thus the isospin transport in heavy-ion re- 
actions depends on both the long-range and the short-range parts of the isospin-dependent 
in-medium nuclear effective interactions, namely, the symmetry potential and the in-medium 
np scatterings cross sections. The former relates directly to the density dependence of the 
symmetry energy E sym (p). To extract information about E sym (p) has been one of the ma- 
jor goals of heavy-ion reactions induced by neutron-rich nuclei [4,5]. Among the potential 
probes proposed so far, the isospin transport has been found very useful for investigating 
the E sym (p) [6-8]. More specifically, using the isospin and momentum dependent transport 
model IBUU04 [9], a symmetry energy of the form E sym (p) rs 32(p/p ) L1 was extracted 
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recently from the MSU data on isospin transport [8]. This conclusion, however, was drawn 
based on transport model calculations using the experimental free-space NN cross sections. 
In this work, we examine effects of the in-medium NN cross sections on isospin transport in 
heavy-ion collisions at intermediate energies. The in-medium NN cross sections are calcu- 
lated within a scaling model according to the nucleon effective masses consistent with the 
nuclear mean field used in the transport model. It is found that the density dependence 
of symmetry energy extracted from the isospin transport data is altered significantly. The 
possibility of determining both the in-medium NN cross sections and the symmetry energy 
corresponding to the same underlying nuclear effective interactions is also discussed. 

II. MOMENTUM DEPENDENCE OF THE MEAN FIELD AND NUCLEON 
EFFECTIVE MASSES IN NEUTRON-RICH MATTER 

The single nucleon potential is one of the most important inputs to all transport models 
for nuclear reactions. Both the isovector (symmetry potential) and isoscalar parts of this 
potential should be momentum dependent due to the non-locality of strong interactions and 
the Pauli exchange effects in many-fermion systems. In the IBUU04 transport model [9] , we 
use a single nucleon potential derived within the Hartree-Fock approach using a modified 
Gogny effective interaction (MDI) [11], i.e., 



U(p,S, p, t, x) = A u (x) — + A t (x) 




(1) 



4 




FIG. 1. (Color online) Symmetry energy as a function of density for the MDI interaction with 
x = 1,0, —1 and —2. Taken from ref. [8]. 

Here 5 = (p n — p p ) / p is the isospin asymmetry of the nuclear medium. In the 
above r = 1/2 (—1/2) for neutrons (protons) and r ^ r'; a — 4/3; f T (f,p) is the 
phase space distribution function at coordinate r and momentum p. The parameters 
A u (x), Ai(x), B, C TiT , C TtT > and A were obtained by fitting the momentum-dependence of the 
U(p,S,p,T,x) to that predicted by the Gogny Hartree-Fock and/or the Brueckner-Hartree- 
Fock (BHF) calculations [12], the saturation properties of symmetric nuclear matter and the 
symmetry energy of about 30 MeV at normal nuclear matter density po = 0.16 fm~ 3 [11]. 
The incompressibility Kq of symmetric nuclear matter at po is set to be 211 MeV consistent 
with the latest conclusion from studying giant resonances [13-15]. The parameters A u (x) 
and Ai(x) depend on the x parameter according to 

AJx) = -95.98 - x , AAx) = -120.57 + x . (2) 

u\ i cr + 1 (7+1 w 

The parameter x can be adjusted to mimic predictions on the density dependence of sym- 
metry energy E sym (p) by microscopic and/or phenomenological many-body theories. Shown 
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in Fig. 1 is the density dependence of the symmetry energy for a; = —2, —1, and 1. 
The last two terms in Eq. (1) contain the momentum-dependence of the single-particle 
potential. The momentum dependence of the symmetry potential stems from the differ- 
ent interaction strength parameters C T y and C T)T for a nucleon of isospin r interacting, 
respectively, with unlike and like nucleons in the background fields. More specifically, we 
use C un iike = —103.4 MeV and Cute = —11.7 MeV. With these parameters, the isoscalar 
potential estimated from (U neutr on + U pro ton)/^ agrees reasonably well with predictions from 
the variational many-body theory [16], the more advanced BHF approach [17] including 
three-body forces and the Dirac-Brueckner-Hartree-Fock (DBHF) calculations [18] in broad 
ranges of density and momentum. 

What is particularly interesting and important for nuclear reactions induced by neutron- 
rich nuclei is the isovector (symmetry) potential. The strength of this potential can be 
estimated very accurately from {Uneutron — U pro ton)/^ [9]. In Fig. 2 the strength of the 
symmetry potential for four x parameters is displayed as a function of momentum and 
density. Here we have only plotted the symmetry potential at sub-saturation densities 
most relevant to heavy-ion reactions studies at intermediate energies. It is noticed that the 
momentum dependence of the symmetry potential is independent of the parameter x. This 
is because the x appears only in the density dependent part of the single nucleon potential of 
Eq. (1) by construction. Systematic analyses of a large number of nucleon-nucleus scattering 
experiments and (p,n) charge exchange reactions at beam energies below about 100 MeV 
indicate undoubtedly that the symmetry potential at po? i- e -? the Lane potential, decreases 
approximately linearly with increasing beam energy E^ in . The data can be well described 
by using U Lane = a - bE kin where a ~ 22 - 34 MeV and b ~ 0.1 - 0.2 [19-22]. One 
should comment that although the uncertainties in both the parameters a and b are large, 
the decreasing feature of the Lane potential with increasing beam energy is very certain. 
This provides a stringent constraint on the symmetry potential. The potential in Eq. (1) 
at po satisfies this requirement very well as seen in Fig. 2. What is very uncertain but 
most interesting is the momentum-dependent symmetry potential at abnormal densities. 
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The experimental determination of both the density and momentum dependence of the 
symmetry potential is thus desired. Heavy-ion reactions provide a unique tool in terrestrial 
laboratories to explore the symmetry potential at abnormal densities with varying momenta. 




FIG. 2. (Color online) Symmetry potential as a function of momentum and density for MDI 
interactions with x = 1,0, —1 and —2. 

One characteristic feature of the momentum dependence of the symmetry potential is 
the different effective masses for neutrons and protons in isospin asymmetric nuclear matter, 
i.e., 

rK = ( + m L d_U L \ 
m T [ p dp J 
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By definition, the effective mass normally depends on the density and isospin asymmetry 
of the medium as well as the momentum of the nucleon. Conventionally, however, the 
effective mass at the Fermi momentum p T = p/(r) is most frequently used to characterize 
the momentum dependence of the nuclear potential. In our calculations of the in-medium 
NN cross sections, we use the general definition of Eq. (3). With the potential in Eq. (1), 
since the momentum-dependent part of the nuclear potential is independent of the parameter 
x, the nucleon effective masses are independent of the x parameter too. 
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FIG. 3. (Color online) Nucleon effective masses at the respective Fermi surface in asymmetric 
matter as a function of density (upper window) and isospin asymmetry (lower window). 
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Shown in Fig. 3 are the effective masses of neutrons and protons at their respective Fermi 
surfaces as a function of density (upper window) and isospin asymmetry (lower window). It 
is seen that the effective mass of neutrons is higher than that of protons and the splitting 
between them increases with both the density and isospin asymmetry of the medium. We 
notice here that the momentum dependence of the symmetry potential and the associated 
neutron-proton effective mass splitting is still highly controversial within different approaches 
and/or using different nuclear effective interactions [23-25]. Being phenomenological and 
non-relativistic in nature the neutron-proton effective mass splitting in the present study is 
consistent with predictions of all non-relativistic microscopic models, see, e.g., [12,17,26], and 
the non-relativistic limit of microscopic relativistic many-body theories, see, e.g., [18,27,28]. 
Recent transport model studies indicate that the neutron/proton ratio at high transverse 
momenta and/or rapidities is a potentially useful probe of the neutron-proton effective mass 
splitting in neutron-rich matter [9,29]. In this work we explore effects of the neutron-proton 
effective mass splitting on in-medium NN cross sections in neutron-rich matter. Applications 
of these in-medium NN cross sections in heavy-ion reactions may be useful for finding other 
probes of the nucleon effective masses in neutron-rich matter. 



III. NUCLEON-NUCLEON CROSS SECTIONS IN NEUTRON-RICH MATTER 

While much attention has been given to finding experimental observables constraining 
the symmetry energy, little effort has been made so far to study the NN cross sections in 
isospin asymmetric nuclear matter. Most of the existing work on the in-medium NN cross 
sections have concentrated on their density dependence in isospin symmetric nuclear matter, 
see, e.g., [30-39]. Here we extend the effective mass scaling model for in-medium NN cross 
sections [30,31,34] to isospin asymmetric matter. Both the incoming current in the initial 
state and the level density of the final state in NN scatterings depend on the effective masses 
of the colliding nucleons. Assuming all matrix elements of the NN interactions are the same 
in free-space and in the medium, the NN cross sections in the medium a^ lum are expected 
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to be reduced compared with their free-space values crj^f by a factor 

Rmedium = &NN = (t^NN I ' ^Nn) i (4) 

where the unn and fi* NN are the reduced masses of the colliding nucleon pairs in free-space 
and in the medium, respectively. 

At relative momenta less than about 240 MeV/c and densities less than about 2p , the 
scaling of a% e / um / 

a NN m Eq. (4) was recently found to be consistent with calculations 
based on the DBHF theory [40]. This finding lends a strong support to the effective mass 
scaling model of the in-medium NN cross sections in the limited density and momentum 
ranges. In this work, we apply the scaling model to elastic NN scatterings in heavy-ion 
reactions at beam energies up to about the pion production threshold. At higher energies, 
inelastic reaction channels become important and in-medium effects on these channels have 
been a subject of much interest, see, e.g., refs. [41-43]. For the inelastic channels we keep 
using the experimental free-space cross sections. This assumption has no effect on our 
present study mainly at intermediate energies. 

While the effective masses and the in-medium NN cross sections have to be calculated 
dynamically in the evolving environment created during heavy-ion reactions, it is instructive 
to examine the in-medium NN cross sections in isospin asymmetric nuclear matter at zero 
temperature. In this situation the integrals in Eq. (1) can be analytically carried out. More 
specifically [11,44], 



,3 / U{r,p) = A A 3 
P l + (p-p>) 2 /A 2 h 3 * 



2pA n (p-Pf(r)) 2 + A 2 



+^ - 2{arctan^M - arctan^^} 
A A A 



(5) 

The medium reduction factor R me dium can thus also be obtained analytically, albeit lengthy. 
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FIG. 4. (Color online) The reduction factor of the in-medium nucleon-nucleon cross sections 
compared to their free-space values as a function of density (top), isospin asymmetry (middle) and 
momentum (bottom). 

As an illustration of a simplified case, we show in Fig. 4 the reduction factor R me dium for 
two colliding nucleons having the same momentum p. The R m edium factor is examined as a 
function of density (upper), isospin asymmetry (middle) and the momentum (bottom). It is 
interesting to note that the in-medium NN cross sections are not only reduced compared with 
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their free-space values, but the nn and pp cross sections are also split while their free-space 
cross sections are the same. Moreover, the difference between the nn and pp scattering cross 
sections grows in more asymmetric matter. The higher in-medium cross sections for nn than 
for pp are completely due to the positive neutron-proton effective mass splitting with the 
effective interaction used. This feature may serve as a probe of the neutron-proton effective 
mass splitting in neutron-rich matter. This possibility will be studied in a future work. We 
also note that the in-medium NN cross sections are also independent of the parameter x. 
They are solely determined by the momentum dependence of the nuclear potential used in 
the model. 



IV. DYNAMICAL GENERATION OF NUCLEON EFFECTIVE MASSES 
DURING HEAVY-ION REACTIONS 

The nucleon effective masses change dynamically in heavy-ion collisions. How big are the 
nucleon effective masses and how much the NN cross sections are modified compared with 
their free-space values in a typical heavy-ion reaction at intermediate energies? To answer 
these questions, we show in Fig. 5 the correlation between the average nucleon effective 
mass and the average nucleon density (top), and the distribution of nucleon effective masses 
(bottom) at the instant of 10 fm/c in the reaction of 132 Sn+ 124 Sn at a beam energy of 50 
MeV/A and an impact parameter of 5 fm. For this particular calculation x = is used, 
very similar results are obtained using other values for the x parameter. It is seen that the 
nucleon effective masses decrease with increasing density. The maximum density reached 
at the instant considered, i.e., 10 fm/c, is about lAp/p . Moreover, the neutron-proton 
effective mass splitting is seen to increase slightly at supra-normal densities. However, the 
increase is not much because the isospin asymmetry normally decreases with increasing 
density, i.e., the isospin fractionation (distillation), see, e.g., [4,5,55]. These features are 
consistent with our expectations discussed in previous sections. From the lower window it 
is seen that the distribution of nucleon effective masses picks at about 0.7 GeV. It is also 
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noticed that some small number of nucleons obtain effective masses above their free masses. 
This is also understandable. In principle, the slope of the nucleon potential du/dp in Eq. 
(3) can be negative during heavy-ion reactions although it happens very rarely, leading to 
the higher effective masses of some nucleons. 




m* T (GeV) 



FIG. 5. (Color online) The correlation between the average nucleon effective mass and the 
average nucleon density (top), and the distribution of nucleon effective masses (bottom) in the 
reaction of 132 Sn+ 124 Sn at a beam energy of 50 MeV/A and an impact parameter of 5fm. 
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FIG. 6. (Color online) The distribution of the reduction factor of in-medium NN cross sections 
in the reaction of 132 Sn+ 124 Sn at a beam energy of 50 MeV/A and an impact parameter of 5fm at 
10, 50 and 100 fm/c, respectively. The insert is the evolution of the central density in the reaction. 

With the nucleon effective masses available we can now examine their effects on nucleon- 
nucleon scatterings during heavy-ion reactions. Shown in Fig. 6 are the distributions of the 
reduction factor R me dium in the reaction of 132 Sn+ 124 Sn at a beam energy of 50 MeV/A and 
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an impact parameter of 5 fm at 10, 50 and 100 fm/c, respectively. The insert in the bottom 
window shows the evolution of the central density during the reaction. The three instants 
represent the compression, expansion and freeze-out stages of the reaction. The quantity 
N pair (Ar < 2.5/m) is the number of nucleon pairs with spatial separations less than 2.5 
fm. These are potential colliding nucleons whose scattering cross section will be reduced by 
the factor R me dium, i-e., <7™^f mm = R me dium x &nn- ^ i s seen that on average as much as 
50% reduction occurs for NN scatterings in the early stage of the reaction. As the system 
expands the average density decreases, the reduction factor Rmedium thus gradually shifts 
towards 1 in the later stage of the reaction. 

V. EFFECTS OF THE IN-MEDIUM NN CROSS SECTIONS ON ISOSPIN 
TRANSPORT IN HEAVY-ION REACTIONS AT INTERMEDIATE ENERGIES 

The in-medium NN cross sections are expected to affect several aspects of heavy-ion 
reactions. Previous studies have already found that the in-medium NN cross sections in- 
fluence both the degree and rate of isospin equilibrium [45-50], see, e.g., [51] for a review. 
However, none of them is in the context of extracting the symmetry energy. In this section 
we investigate how the symmetry energy extracted from the isospin transport data might 
be altered. 



A. Isospin diffusion/transport 

Experimentally, one way of measuring quantitatively the amount of isospin transport, 
between the projectile nucleus A and the target nucleus B is to study the quantity Ri defined 
as [52] 

2X A+B - X A+A - X B+B 

^ ~ x A+A — X B+B 

where X is any isospin-sensitive observable. The Ri is also known as a useful tool of measur- 
ing quantitatively the degree of isospin diffusion in heavy-ion reactions [7] . Here we prefer to 
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use the more general term isospin transport although the two terminologies may have been 
used interchangeably in the literature. Our main concern is that the term diffusion is often 
associated with irreversible processes. While by definition as given in Eq. (6), and from 
analyzing its time evolution in the following, it is seen that the Ri does not always have to 
decrease with time. In other words, there is no guarantee that the particular observable X is 
synchronized in all three reaction systems such that Ri always decreases with time. This is 
because the evolution and freeze-out time of the observable X may depend significantly on 
the system size especially if the masses of A and B are rather different. By construction, the 
value of Ri is 1 (—1) for the symmetric A + A (B + B) reaction. If an isospin equilibrium is 
reached as a result of isospin transport the value of Ri is about zero. As mentioned earlier, 
the EOS of symmetric nuclear matter still has some uncertainties associated with the density 
and momentum dependence of the isoscalar potential. Fortunately, these uncertainties and 
all others due to the isospin-independent ingredients in the reaction dynamics can be largely 
canceled out because of the special construction of Ri. The Ri also has the advantage of 
minimizing significantly effects of pre-equilibrium emissions [7]. 
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FIG. 7. (Color online) The evolution of isospin diffusion Ri using the four x parameters and 
the in- medium nucleon-nucleon cross sections. 



16 



In the NSCL/MSU experiments with A = 124 Sn on B = 112 Sn at a beam energy of 50 
MeV/nucleon and an impact parameter about 6 fm, the isospin asymmetry of the projectile- 
like residue was used as the isospin tracer X [7]. The data was recently analyzed within the 
IBUU04 transport model using the free-space experimental NN cross sections. Consistent 
with the experimental selection, in model analyses the average isospin asymmetry (5) of the 
124 Sn-like residue was calculated from nucleons with local densities higher than po/20 and 
velocities larger than 1/2 the beam velocity in the cm. frame. Shown in Fig. 7 are the 
time evolutions of Ri re-calculated using the in-medium NN cross sections and the four x 
parameters. The data from MSU is indicated by the shaded band. It is seen that the net 
isospin transport and the influence of the x parameter show up mainly in the expansion 
phase of the reaction after about 40 fm/c. The values of Ri stabilize approximately after 
about 80 fm/c. It is seen that the values of Ri with x — — 1 and come close to the MSU 
data in the late stage of the reaction. 

For a more meaningful comparison with the experimental data, we have calculated the 
time average of Ri between t = 120 fm/c and 150 fm/c as in Ref. [8]. Shown in Fig. 8 is 
a comparison of the averaged strength of isospin transport 1 — Ri obtained with the free 
and in-medium NN cross sections, respectively, as a function of the asymmetric part of the 
isobaric incompressibility of nuclear matter at po [53-55] 

K asy ( Po ) = 9p 2 (d 2 E sym /dp 2 )^ - 18p (dE sym /dp) p() . (7) 

In each case, 2000 events were generated for all three reaction systems used in the analysis. 
The error bars were drawn to indicate fluctuations and were obtained from the dispersion 
of Ri time evolution [8]. First, it is interesting to note that with the in-medium NN cross 
sections the strength of isospin transport 1 — R^ decreases monotonically with the decreasing 
value of x. With the free-space NN cross sections, however, there appears to be a minimum 
at around x — — 1. Moreover, this minimum is the point closest to the experimental data. 
This allowed us to extract the value of K asy (p ) = —550 ± 100 MeV. With the in- medium 
NN cross sections we can now further narrow down the K asy (po) to be about —500 ± 50 
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MeV. The latter is consistent with that extracted from studying the isospin dependence of 
giant resonances of 112 Sn to 124 Sn isotopes by Fujiwara et al at Osaka [56]. Shown also in 
the figure are the 7 values used in fitting the symmetry energy with E sym (p) = 31.6(p/po) 7 - 
The results with the in-medium NN cross sections constrain the 7 parameter to be between 
0.69 and 1.05. The lower value is close to what is extracted from studying giant resonances 
[57,58]. The value of 7 = 1.05 extracted earlier using the free-space NN cross sections sets 
an upper limit. 
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FIG. 8. (Color online) The degree of isospin transport as a function of K asy {po) with the free 
(filled squares) and in-medium (open squares) nucleon-nucleon cross sections. 



It is seen that the difference in 1 — Ri obtained with the free-space and the in-medium NN 
cross sections is about the same with x = 1 and x = 0, but then becomes especially large at 
x = — 1 and x = —2. Why does the effect of the in-medium NN cross sections increase with 
the decreasing K asy (p ) or x parameter? This question can be understood from considering 
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contributions from the symmetry potential and the np scatterings. As we have mentioned 
in the introduction, both contributions to the isospin transport depend on the np scattering 
cross section o np . Schematically, the mean field contribution is proportional to the product 
of the isospin asymmetric force F np and the inverse of the np scattering cross section a np . 
While the collisional contribution is proportional to the a np . The overall effect of the in- 
medium NN cross sections on isospin transport is a result of a complicated combination of 
both the mean field and the NN scatterings. Generally speaking, the symmetry potential 
effects on the isospin transport will become weaker when the NN cross sections are larger 
while the symmetry potential effects will show up more clearly if smaller NN cross sections 
are used. This feature can be seen from Fig. 9. 
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FIG. 9. (Color online) Symmetry potential as a function of momentum at selected densities 
(a), and density at selected momenta (b) with the four x parameters. Taken from Ref. [8]. 

For x = 1 and x = 0, the symmetry potential and its density slope, as shown in Fig. 9, 
are large at low densities where the majority of net isospin transport occurs. The F np factor 
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makes the contribution due to the mean field dominates over that due to the collisions. 
Therefore, the reduced in-medium a np leads to about the same or a slightly higher isospin 
transport. As the x parameter decreases to x — — 1 and x = —2, however, the symmetry 
potential decreases and its density slope can be even negative at low densities. Thus in these 
cases either the collisional contribution dominates or the mean field contribution becomes 
negative. The reduced in-medium np scattering cross section a np leads then to a lower 
isospin transport compared with the case with the free-space NN cross sections. 
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FIG. 10. (Color online) The skin thickness in lead, 5R, versus the isospin diffusion parameter 
1 — Ri for the four equations of state in this work. The present acceptable experimental range for 
5R and the NSCL/MSU data on isospin diffusion are also shown. 



20 



B. Correlation between the isospin diffusion and neutron-skin in 08 Pb 



Both the neutron skin of heavy nuclei and the degree of isospin diffusion in heavy-ion 
collisions at intermediate energies are sensitive to the symmetry energy at sub-saturation 
densities. Since they both are determined by the same underlying EOS, the correlation 
between these seemingly very different observables provides a more stringent constraint on 
the symmetry energy than their individual values [59] . In view of the revised calculations of 
the isospin diffusion using the in-medium NN cross sections, here we revisit the correlation 
between the isospin diffusion and the size of neutron-skin in 208 Pb. Shown in Fig. 10 is 
the isospin diffusion parameter I — Ri versus the skin thickness in lead, SR. The presently 
acceptable experimental range for 5R [60] and the NSCL/MSU data on isospin diffusion are 
also shown. The 5R, taken from Ref . [59] , was calculated within the Hartree-Fock approach 
using the same equations of state corresponding to the four x parameters. While the isospin 
diffusion data is between calculations with x = and x — — 1, the upper limit of the present 
n-skin measurements comes closer to the prediction using x — 0. On one hand, the present 
analysis of the isospin diffusion data favors a neutron skin as large as about 0.25 fm. On 
the other hand, the present neutron-skin data favors x = and x — 1. One can thus 
conclude that a symmetry energy of E sym (p) ps 31.6(p/p ) ' 69 is currently most acceptable 
based on both the neutron-skin and the isospin diffusion data. It is seen that more accurate 
measurements of both kinds of experimental data, especially the neutron-skin, are very 
desirable. 



VI. NUCLEON TRANSVERSE FLOW AS A PROBE OF THE IN-MEDIUM NN 

CROSS SECTIONS 

Our results above indicate clearly that the in-medium NN cross sections affect signif- 
icantly the extraction of symmetry energy from isospin transport in heavy-ion reactions. 
Therefore, auxiliary measurements of other observables have to be made to constrain the 
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in-medium NN cross sections, preferably from the same experiments. Fortunately, a number 
of observables are known to be sensitive to the in-medium NN cross sections. These in- 
clude the quadruple moment Q zz of nucleon momentum distribution, the linear momentum 
transfer (LMT) and the ratio of transverse to longitudinal energies (ERAT), see, e.g., ref. 
[39,43,61,62]. It is also well known that these observables are also sensitive to the EOS of 
symmetric matter. Therefore, the extraction of the in-medium NN cross sections from these 
observables relies closely on our knowledge about the EOS of symmetric nuclear matter. 
Nevertheless, great progress has been made over the last three decades in determining the 
EOS of symmetric nuclear matter. Using a combination of several observables in heavy- ion 
reactions, such as the elliptic flow and kaon production, see, e.g., Refs. [63-65], and the 
analysis of giant resonances, the EOS of symmetric nuclear matter has been severely con- 
strained. It has now been widely recognized that the major remaining uncertainty in further 
constraining the EOS of symmetric nuclear matter is our poor knowledge about the symme- 
try energy [14,15,64]. In this work we thus concentrate on the possibility of extracting the 
symmetry energy and the in-medium NN cross sections simultaneously without considering 
the remaining uncertainty in the EOS of symmetric matter. 

In this work we choose to use the transverse flow to constrain the in-medium NN cross 
sections. Its sensitivity to variations of the in-medium NN cross sections is well known, 
especially around the balance energies, see, e.g., [66-73]. As an illustration, shown in Fig. 
11 is the transverse flow of all free nucleons identified as those having local densities less 
than po/8 & t freeze-out in the reaction of 132 Sn+ 124 Sn at an impact parameter of 5 fm and 
a beam energy of 400 MeV/A (upper) and 50 MeV/A (bottom) with the free-space and 
in-medium NN cross sections, respectively. It is seen that at 50 MeV/A the transverse flow 
is much more sensitive to the in-medium NN cross sections than at 400 MeV/A. In fact, 
the direction of transverse flow is altered from being positive to slightly negative by the 
reduced in-medium NN cross sections. This observation is consistent with our expectation 
and results of previous studies. 



22 



250 



132 Sn+ 124 Sn, b=5 fm 




CD 



e- 50 , , , , | , , , , | , , , , | , , , , | , , , , | , , , , 

A 

< 40 - 

X 

v 30 - 




-40 - 

-50 I 

-1.5 -1 -0.5 0.5 1 1.5 



cms j beam 



FIG. 11. Nucleoli transverse flow in the reaction of Sn+ Sn at an impact parameter of 5 
fm and a beam energy of 400 MeV/A (upper) and 50 MeV/A (bottom). 
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FIG. 12. Neutron transverse flow in the reaction of 112 Sn+ 112 Sn (left), 124 Sn+ 124 Sn (middle) 
and 124 Sn + 112 Sn (right) at an impact parameter of 6 fm and a beam energy of 50 MeV/A. 
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FIG. 13. Proton transverse flow in the reaction of 112 Sn + 112 Sn (left), 124 Sn + 124 Sn (middle) 
and 124 SVi + 112 Sn (right) at an impact parameter of 6 fm and a beam energy of 50 MeV/A. 
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Shown in Figs. 12 and 13 are the transverse flow for all (free and bound) neutrons and 
protons, respectively, with x = and x = —1 for the three reactions used in studying 
the isospin transport. Free nucleons show the same features but with poor statistics for 
the same number of events. There is very little system size dependence among the three 
reactions considered. It is seen that the transverse flow is much more sensitive to the in- 
medium NN cross sections than to the symmetry energy parameter x. There is a weak 
sensitivity to the variation of the x parameter, especially with the reduced in-medium NN 
cross sections for neutrons. Protons are affected by both the repulsive Coulomb potential 
and the generally attractive symmetry potential, while neutrons are only affected by the 
repulsive symmetry potential, besides the same isoscalar potential acting on both neutrons 
and protons. Neutrons thus appear to be more sensitive to the x parameter. Moreover, the 
Coulomb potential dominates over the symmetry potential, leading to the slightly higher 
transverse flow for protons than for neutrons as seen from comparing Fig. 12 with Fig. 
13. It is also interesting to note from comparing Fig. 12 with Fig. 13 that there is a clear 
indication of higher (lower) transverse flow for neutrons (protons) with x = than that 
with x = —1, especially around the projectile and target rapidities. This feature is what 
one expects from considering the symmetry potentials shown in Fig. 9. At sub-saturation 
densities the symmetry potential is stronger with x = than that with x — — 1. Therefore, 
the neutron (proton) transverse flow is stronger (weaker) with x — 0. This fine dependence 
on the symmetry energy parameter x can be studied in more detail by using the neutron- 
proton differential flow [74]. We find, however, that the strength of the neutron-proton 
differential flow only changes by about 1 to 2 MeV/c by varying the x parameter from —1 to 
0. It is much less than the change due to the variation of the in-medium NN cross sections. 
Based on these results, we therefore propose to measure the proton transverse flow in order 
to constrain the in-medium NN cross sections. 
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FIG. 14. Combined nucleon transverse flow in the reaction of 124 Sn + 124 Sn at an impact 
parameter of 6 fm and a beam energy of 50 MeV/ A. 

To go one step further, taking into account the generally positive/negative nature of 
the symmetry potentials for neutrons/protons, one may use the combined transverse flow 
of neutrons and protons < p^ p /A > to better determine the in- medium NN cross sections 
without much influence of the symmetry potential. This, of course, requires measuring 
neutrons as accurately as protons simultaneously. Shown in Fig. 14, is an example of the 
combined nucleon flow for the 124 Sn+ 124 Sn reactions. It is clearly seen that the combined 
flow of all nucleons are much less affected by the variation of the symmetry energy, making 
it an even better probe of the in-medium NN cross sections. 



VII. SUMMARY 

In order to better understand the isospin dependence of the in-medium nuclear effective 
interactions, we investigated effects of the in-medium NN cross sections on isospin trans- 
port in heavy-ion reactions within the transport model IBUU04. The isospin-dependent 
in-medium NN cross sections consistent with the nuclear mean field used in the transport 
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model were evaluated by using the scaling model according to the nucleon effective masses. 
It is found that the NN cross sections in neutron-rich matter are not only reduced compared 
with their values in free space, their isospin dependence is also altered. Because of the posi- 
tive and growing neutron-proton effective mass splitting in more neutron-rich matter for the 
effective interactions used in this work, the splitting between the nn and pp cross sections 
increases with the increasing isospin asymmetry of the medium. 

The in-medium NN cross sections are found to influence significantly the isospin transport 
and nucleon transverse flow in heavy-ion reactions at intermediate energies. By using the 
free-space experimental NN cross sections, a symmetry energy of E sym (p) pa 31.6(p/po) 1 ' 1 
was extracted from the MSU data on isospin transport. With the in- medium NN cross 
sections, however, the symmetry energy of E sym (p) pa 31.6(p/p ) ' 69 was found most accept- 
able in comparison with both the MSU isospin diffusion data and the presently acceptable 
neutron-skin thickness in 208 Pb. The isospin dependent part K asy (po) of the isobaric im- 
compressibility of nuclear matter was further narrowed down to —500 ± 50 MeV. 

The possibility of determining simultaneously both the in-medium NN cross sections and 
the symmetry energy corresponding to the same underlying nuclear effective interactions was 
also studied. The proton transverse flow, or even better the combined transverse flow of 
neutrons and protons all together, can be used as an effective probe of the in-medium NN 
cross sections without much hindrance from the uncertainties of the symmetry energy. Our 
findings in this work demonstrated clearly the importance of using the in-medium NN cross 
sections consistent with the momentum-dependent nuclear mean field in transport model 
studies of heavy-ion reactions. 
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